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Recent results on moments of parton distribution functions * 
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We report on recent results for the pion and nucleon matrix element of the twist-2 operator corresponding to 
the average momentum of non-singlet quark densities. We discuss finite size effects for the nucleon matrix element 
and present first preliminary results for the non-perturbative renormalisation from full dynamical simulations. 



1. INTRODUCTION 

Phcnomenological fits to experimental data 
provide results for moments of parton distribu- 
tion functions which allow a direct comparison 
with lattice calculations. In a recent paper P 
we were able to compute a continuum limit value 
of the lowest moment of a twist-2 operator, corre- 
sponding to the average momentum of non-singlet 
quark densities, in pion states. In obtaining the 
renormalisation group invariant matrix element, 
we have controlled important systematic errors, 
such as non-perturbative renormalisation [2], fi- 
nite size effects [H] and effects of a non- vanishing 
lattice spacing. The crucial limitation of our cal- 
culation was the use of the quenched approxima- 
tion. 

In the scope of our current investigations is 
the comparison of finite size effects for pion and 
nucleon matrix elements. In addition we con- 
centrate on computing non-perturbatively the 
evolution of the twist-2 operators, using the 
Schrodinger functional finite-size technique, on 
Nf = 2 dynamical configurations provided by the 
ALPHA collaboration H]. 

2. FINITE SIZE EFFECTS 

Finite size effects for nucleon matrix elements 
have not been studied in detail before, but might 
infiuence considerably the precise determination 
of moments of parton distribution functions from 
lattice QCD. We have investigated the effects 
of limited lattice extent at (3=6.0 (a ~ 0.093 
fm) with the non-perturbatively improved clover 
fermion action and lattice sizes (12 — 32)^ x 32. 

The twist-2 non-singlet operator for a flavour 



doublet of fermions belongs to two irreducible 
representations of the lattice H{A) group 

Ci2(x) = ^Vi(x)7{i D2} T^V'(x) 
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For the study of finite size effects (FSE) we con- 
centrate solely on the ©44 operator, since it can 
be computed at zero external momentum and 
thus provides a reliable signal. In particular we 
use Schrodinger functional (SF) boundary condi- 
tions (0J for our computation, which allows the 
extraction of the matrix element from one large 
plateau around the middle of the time extent T of 
the lattice. The correlation function of the matrix 
element is obtained by inserting the 044 op- 
erator between two SF nucleon states § and S' at 
the time boundaries t = Q and T and suitable nor- 
malisation with the boundary-to-boundary corre- 
lator Ci 
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The bare nucleon matrix element is defined by 
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= {N\Oi4N) = .fM{xo/L)/C, . 
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In order to obtain the infinite volume limit 
nucleon mass and matrix element we performed 
a purely phenomenological fit X{L) = cq + 
ci/L'^/^ exp(— C2i), where X^uifqa or {x)u-d- 
Please note that a power-law fit ansatz would de- 
scribe the data almost equally well. In figs.^and 
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Figure 1. Finite size effects for m^a 



R<x>„(L) = (<x>N(L)-<x>fj)/<x>N 
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Figure 2. Finite size effects for nucleon {x)u-d 

El we show tlie relative systematic deviation from 
the infinite volume limit values for both the nu- 
cleon mass and matrix element. While the nu- 
cleon mass shows a deviation of about 4% at 
m^L ~ 5, the nucleon matrix element is over- 
estimated by 10-30% at this value of m^i. In 
both cases one can be sure to exclude FSE only 
for tjIt^L > 7. For the pion matrix element we ob- 
served a FSE of about 3% at tjIt^L ~ 5, while the 
pion mass is found to be free from FSE already 

for TTlTrL > 4. 

These observations lead to the expectation that 
the strong FSE for nucleon matrix elements could 
be one ingredient for the deviation of current lat- 
tice QCD estimates of {x)u-d in comparison with 
experimental data. 



3. Nf = 2 CONTINUUM STEP SCALING 
FUNCTION 

The strategy used to compute the non- 
perturbative evolution of the operators ©44 and 
has been described in detail in our study 
for the quenched case The evolution from 
initially large L/low /i to small L/high /i is ob- 
tained by applying the so called step scaling func- 
tion (SSF) in the Schrodinger functional scheme. 
Once the perturbative regime is safely reached 
one continues the evolution in perturbation the- 
ory computing the scale and scheme independent 
RGI matrix element. The connection with exper- 
imental results can be obtained by making the 
adequate perturbative evolution of the RGI ma- 
trix element in the MS scheme. 

The renormalisation conditions for the local 
044 and O12 operators are given by 

Onifi) = Z^\afi)0{a), Onifi = 1'^) = . 

The correlators to compute the Z factor are 

fo{x,^/L) = -\ J2 (0(x)§,(y,z)) 



fi^-l E (S;(u,v)S,(y,z)) 

u.v.y .z 

where Sq and are suitable quark sources to 
probe the operators O. With this definition the 
renormalisation constants are obtained by 



Z[a/L,^)~c p= — ; c 



In order to map the L dependence recursively we 
use the SSF, rigorously defined on the lattice by 



lim Y^Zo (iJ-^ o/L) 



Zoju, 2L/a) 
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The values of (3 corresponding to a fixed running 
coupling for Nf = 2 dynamical 0(a)-improved 
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Figure 3. Continuum step scaling function (7z4,4, 

Wilson fermions are available in |4l5j . We have 
computed the SSF at six values of the renor- 
malised couphng g%p{L) = 0.98 to 3.33 for the 
lattice sizes L/a = 6,8,12, which includes the 
determination of Z factors on the corresponding 
L/a and 2L/a lattices. The continuum limit has 
been obtained by a weighted average of several 
independent runs on the larger L/a = 8, 12 (and 
2L/a) lattices using the HMC with one or two 
pseudo-fermion fields j7j . In fig. O we show the 
A^^ — 2 continuum running of cr^^^ in comparison 
with the previously obtained quenched evolution 
. In addition we plot the universal 1-loop (dot- 
ted) and 2-loop order {Nf — 0, dashed) results 
of perturbation theory. The corresponding plot 
for the running of (yzi2 shows a very similar be- 
haviour. 

Complemented by a computation of the pion 
or nuclcon matrix element for Nf = 2 dynamical 
fermions, the SSF allows to obtain the RGI ma- 
trix element completely non-perturbatively. The 
whole strategy is summarised by the following for- 
mula: 
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where we use the n — 6 SSF computed at /i„ 

a(^,g^(/.))=.(^,g^(Mi))---a(^,g^(M„)) 
Mo Mo Mn-l 

to jump from the non-perturbative scale fiQ to the 



perturbative (ultraviolet) scale fi. At this point 
one can try to do the perturbative matching using 
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computed with 3-loop f3 and 2-loop 7 functions. 
If the perturbative matching has been successful 
the quantity 

should be independent from the ultraviolet scale 
M- Indeed we found in the quenched case that the 
last steps give a very nice plateau ■ The deter- 
mination of cT/'/vV for Nf = 2 dynamical fermions 
requires a new computation of the 2-loop anoma- 
lous dimension in the Schrodinger functional 
scheme, which is currently work in progress. 
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